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ABSTRACT 

A  procedure  has  been  developed  for  processing  the  strains 
available  within  the  computer  program  ADINA  to  provide  strain 
output  in  any  desired  orientation  and. to  calculate  and  output 
strains  on  the  surface  of  a  three-dimensional  solid  or  shell 
element.  To  this  end,  prpper  additions  and  modifications  have 
been  made  to  ADINA  to  implUnfctit  the  procedure  directly  within 
the  program  to  allow  for  thp  processing  and  output  of  strains 
during  the  normal  execution  stage.  This  includes  reading  in 
additional  input  data  (of  minimum  Amount),  transformation  of 
strains  from  the  global  system  to  specified  coordinate  direc¬ 
tions,  the  output  of  transformed  strains,  and  the  calculation 
and  output  of  surface  strains  based  on  the  interpolation 
functions  available  in  ADINA.  At  the  present  time  such  capa¬ 
bilities  are  provided  for  a  limited  number  of  elements  on  a 
selective  basis,  in  order  to  minimize  the  amount  of  computer 
storage  required. 

Some  typical  structures  have  been  analyzed  to  compute  and 
output  strains  both  at  interior  points  ana  on  exterior  surfaces 
of  elements,  for  a  set  of  elements  specified  by  the  user  during 
the  input  stage.  Results  demonstrate  the  feasibility  of  the 
procedure  fo_-  the  transformation  and  extrapolation  of  strains. 

ADMINISTRATIVE  INFORMATION 

The  work  reported  herein  was  funded  by  the  Naval  Sea  Systems  Command  (NAVSEA 
032R),  under  the  Exploratory  Development  Program  in  Submarine  Combat  Protection.  The 
work  was  performed  at  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
under  Program  Element  62543N,  Task  Area  SF43400001,  and  Work  Unit  1750-294. 

Dr.  Theodore  G.  Toridis  is  a  Professor,  School  of  Engineering  and  Applied 
Science,  George  Washington  University,  Washington,  D.C.  20052. 

INTRODUCTION 

The  analytical  solution  of  most  practical  structures  culminates  in  the  computa¬ 
tion  of  displacements,  stresses  and  strains,  or  other  quantities  of  interest  to  the 
engineer.  Whenever  efforts  are  made  to  compare  analytical  results  with  test 
measurements,  the  comparison  of  analytical  and  experimental  strains  plays  an  impor¬ 
tant  role  in  the  process  of  verification  of  a  structural  model,  because,  in  general, 
strain  measurements  are  relatively  easy  to  obtain  and  record  during  the  testing 
process . 


When  a  material  behaves  elastically,  strain,  in  the  engineering  sense,  is 
ordinarily  defined  as  the  change  in  length  per  uhit  initial  length.  This  definition 
of  strain  was  first  introduced  by  Cauchy  and  is  usually  measured  quantitatively  by 
performing  a  uniaxial  tension  test  and  determining  the  changes  in  length  as  the 
specimen  deforms.  When  finite  strains  and  deformations  take  place  in  a  body  in  a 
three-dimensional  sense,  the  problem  of  characterising  strain  becomes  more  compli¬ 
cated.  In  such  cases,  two  different  definitions  of  strain  are  commonly  employed,  one 
due  t.o  Green  and  the  other  due  to  Almansi.  If  the  strains  are  referred  to  the 
initial  (usually  undeformed)  configuration  of  the  system,  they  correspond  to  the 
so-called  Lagrangian  formulation.  If,  on  the  other  hand,  strains  are  expressed  as 
functions  of  the  deformed  configuration,  they  are  normally  identified  with  the 
Eulerian  formulation.  If  large  deformations  are  associated  with  material  and 
geometric  nonlinearities,  coupled  with  large  strains,  the  computation  of  strains  in 
a  complex  structural  system  becomes  a  formidable  task,  and  can  be  accomplished 
effectively  only  by  the  use  of  a  general  purpose  finite  element  computer  program, 
such  as  ADINA, 

In  ADINA,  in  the  case  of  a  linear  elastic  continuum,  components  of  the  strain 

tensor  t+'^te  at  time  t  +  At  are  computed  incrementally  based  on  the  usual  infini- 
rs 

tesimal  theory  and  expressed  as 
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where  t+AtQur  s  are  the  displacement  gradients  at  time  t+At  measured  with  reference 
to  the  initial  configuration  of  the  body  at  time  0.  In  the  case  of  large  deforma¬ 
tions,  when  the  Total  Lagrangian  Formulation  coupled  with  an  incremental  procedure 
is  used,  the  strain  relations  are  represented  by 
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*A  complete  listing  of  references  is  given  on  page  21. 
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When  the  updated  Lagrangian  Formulation  is  employed,  the  corresponding  expres¬ 
sions  for  the  strains  are  given  as 
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Relations  of  this  type  can  be  employed  to  compute  strains  as  part  of  the  overall 
computational  process.  In  the  case  of  a  thin  shell  element,  complete  strain- 
displacement  transformation  matrices  for  the  ADINA  shell  element,  for  both  the  linear 
and  nonlinear  cases,  have  been  derived  and  given  in  Reference  2. 

The  basic  computational  algorithm  embedded  in  ADINA  makes  provision  for  the 
computation  of  strains,  but  no  special  effort  is  made  to  generate  elaborate  strain 
output  at  the  option  of  the  user  of  the  program.  More  specifically,  the  strain 
output  provided  in  ADINA  allows  for  storing  the  strains  in  a  porthole  for  subsequent 
access,  but  the  strains  are  referred  to  the  global  Cartesian  coordinate  system. 

For  a  linear  material  model,  the  strains  are  available  at  certain  nodes  defined  by 
the  stress  tables,  or  at  the  element  integration  points.  For  nonlinear  material 
models,  the  strains  are  available  at  the  integration  points  only. 
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In  the  present  study,  for  a  three-dimensional  solid  and  shell  element,  the  • 
above  capabilities  were  extended  in  order  to  provide  for  a  more  comprehensive  and 
versatile  strain  output.  At  present,  the  options  provided  are  summarized  belows 

1.  For  linear  material  models,  strains  are  output  at  nodal  points  specified 
by  stress  tables.  However,  the  strains  can  be  referred  to  an  arbitrarily  oriented 
Cartesian  coordinate  system  by  inputting  the  corresponding  direction  cosine  matrix 
for  each  element.  Alternatively,  the  user  can  specify  the  output  of  strains  with 
reference  to  the  natural  (r,s,t)  coordinate  system  as  defined  by  ADINA. 

2.  For  linear  or  nonlinear  material  models,  the  strains  are  output  at 
integration  points  with  reference  to  any  coordinate  system  as  explained  in  Item  1, 
or  to  the  natural  (r,s,t)  coordinate  system. 

3.  For  either  linear  or  nonlinear  material  models,  in  addition  to  the  above 
output,  the  option  of  providing  printout  of  the  strains  on  any  two  surfaces  of  a 
three-dimensional  solid  or  shell  element  has  been  provided.  For  this  purpose,  the 
user  can  specify  a  variable  mesh  for  output  of  strains  at  the  mesh  points  (a  3*3  to 
5x5  mesh),  or  input  an  array  of  coordinate  locations  for  strain  output.  In  computing 
these  strains,  the  interpolation  functions  available  within  ADINA  are  employed.  As 
pointed  out  previously,  this  option  is  significant  in  engineering  research  applica¬ 
tions,  because  experimental  strains  are  measured  on  the  structure  surface  and  must 

be  correlated  with  analytically  computed  strains  at  the  same  locations. 

INTERPOLATION  AND  EXTRAPOLATION  OF  STRAINS 

As  mentioned  earlier,  in  the  case  of  a  linear  material  model  of  a  structure 
made  up  of  three-dimensional  or  shell  elements,  the  output  of  stresses  and  strains 
in  ADINA  is  provided  at  selected  nodal  points  specified  by  means  of  stress  tab3.es. 

The  computation  of  strains  at  any  other  point  in  the  interior  or  at  the  surface  of 
the  element  must  invariably  involve  an  interpolation  process.  Within  the  framework 
of  a  stress  table  employed  in  ADINA,  a  maximum  of  16  points  can  be  employed  to  perform 
such  interpolation.  On  the  other  hand,  for  nonlinear  material  models,  stresses  and 
strains  can  be  output  only  at  Gaussian  integration  points.  Because  the  maximum 
order  of  numerical  integration  in  each  of  the  three  coordinate  directions  is  4,  at 
most  a  total  of  64  points  can  be  used  to  perform  extrapolation. 


Although  such  approximations  can  be  curried  out  based  on  a  variety  of 
procedures,  for  the  purpose  of  accomplishing  this  with  a  moderate  amount  of  effort 
and  a  reasonable  number  of  modifications  to  ADINA,  the  use  of  the  shape  functions 
and  associated  derivatives  inherent  in  ADINA  appears  to  be  the  logicel  approach  to 
follow,  as  was  done  in  the  present  study.  It  should  be  pointed  out,  however,  that, 
in  the  case  of  the  three-dimensional  solid  element  or  shell  element,  the  interpola¬ 
tion  process  involved  in  ADINA  uses  functions  of  order  no  higher  than  two  (quadratic). 
Nevertheless,  in  the  case  of  regular  geometrical  surfaces,  matching  the  location  of 
computed  strains  with  measured  strains  resulting  from  an  experiment  would  ordinarily 
pose  no  difficulty.  On  the  other  hand,  for  bodies  of  arbitrary  shape,  a  more 
accurate  description  of  the  location  of  surface  points  may  be  desirable.  Also,  in 
view  of  the  fact  that  up  to  four  Gaussian  integration  points  can  be  employed  in  each 
coordinate  direction,  the  extrapolation  process  for  the  strains  can  be  based  on 
higher  order  functions. 

3 

In  addition,  as  proposed  by  Hinton  and  Campbell  concerning  an  improved 
extrapolation  method  for  stresses,  in  a  similar  manner  local  least  squares  smoothing 
techniques  can  be  used  to  perform  an  extrapolation  of  strains  from  sampling  locations 
corresponding  to  the  Gaussian  integration  points.  Alternatively,  the  concept  of 

4 

conjugate  approximations  developed  by  Oden  and  Brauchli  for  obtaining  consistent 
approximations  of  stress  fields  in  the  displacement  formulation  of  the  finite  element 
method  can  be  extended  to  the  case  of  approximation  of  strain  fields.  It  should  be 
pointed  out,  however,  that  such  approaches  would  require  major  modifications  to 
ADINA  leading  to  extended  efforts  beyond  the  normal  level, 

Therefore,  from  tne  viewpoint  of  computational  efficiency,  the  use  of  the 
original  shape  functions  employed  in  the  description  of  the  element  geometry  and 
displacements  to  perform  interpolation  or  extrapolation  of  the  strains,  appears  to  be 
the  most  feasible  approach  to  follow  for  achieving  the  immediate  goals  of  this  study. 
Such  a  procedure  can  be  used  to  evaluate  both  linear  and  nonlinear  strains  at  any 
point  in  an  element.  The  distinct  advantage  of  this  approach  is  that  the  process  of 
interpolation  or  extrapolation  for  the  geometry,  the  displacements  and  strains,  are 
all  based  on  the  same  functions  leading  to  a  more  consistent  and  efficient  procedure 
as  far  as  the  computational  effort  is  concerned. 
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USE  OF  ADINA  INTERPOLATION  FUNCTIONS 

The  shape  functions  used  in  ADINA  for  an  8  to  21  vutiable-node-nutnber  three 
dimensional  solid  element  are  in  the  form  of  (see  Figure  1): 
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i  '  ■ 

'■  ,  Similar  shape  functions  are  used  for  the  ADINA  shell  element  except  that  the  more 
geper.al  concept  of  displacement  and  rotation  isoparametric  element  is  used  in  its 
development. ^ 

These  shape  functions  form  the  basis  of  evaluation  of  the  element  properties 
which  are  then  synthesized  to  obtain  the  corresponding  properties  of  the  entire 
structure  taking  into  account  t’.2  boundary  and  constraint  conditions.  Once  the 
displacements  at  the  nodal  points  have  been  solved,  the  computation  of  the  displace¬ 
ments  and,  consequently,  the  strains  anywhere  within  an  element  can  be  accomplished 
using  the  shape  functions  and  associated  derivatives.  As  discussed  previously, 
provisions  were  originally  made  in  ADINA  for  the  computation  of  strains — for  both 
linear  and  nonlinear  material  models — at  specific  points,  with  reference  to  the 
global  coordinate  system.  The  procedure  used  in  computing  these  strains  at  the  nodal 
or  Gaussian  integration  points  has  simply  been  repeated  in  this  study  to  provide 
strain  output  at  other  points  within  the  element  or  on  its  surface.  For  this 
purpose,  it  is  only  necessary  to  specify  the  desired  values  of  the  natural  r,  s,  t 
coordinates  at  a  point  at  which  the  strains  are  to  be  computed.  However,  in  general, 
the  x,  y,  and  z  coordinates  of  the  point  will  be  known,  and  the  determination  of  the 
corresponding  natural  coordinates  for  complicated  geometries  may  involve  an  iterative 
process . 

In  the  case  of  surface  strain  output,  two  different  options  are  provided  as 
follows: 

1.  Surface  strains  are  output  at  points  specified  during  the  input  stage  by 
means  of  an  array  of  natural  (r,s,t)  coordinates  corresponding  to  the  output 
locations . 
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2.  In  order  to  minimize  the  input,  a  variable  mesh  size  is  selected  as  part 
of  the  input,  so  that  surface  strains  can  be  output  automatically  at  locations 
corresponding  to  the  grid  points  of  the  particular  mesh.  To  this  end,  mesh  sizes  of 
3  x  3  to  5  x  5  can  be  used  as  shown  in  Figure  2,  in  which  the  grid  line  locations  are 
also  indicated  by  specifying  the  constant  coordinate  values  along  grid  lines.  These 
locations  have  been  selected  to  correspond  to  Gaussian  integration  points  of 
corresponding  order. 

Once  the  strains  are  computed  at  each  point  of  interest,  they  are  transformed 
from  the  global  coordinate  directions  to  any  desired  direction  using  a  standard 
rotation  matrix  [X],  The  entries  of  [X]  are  the  direction  cosines  associated  with 
the  global  (x^.x^jX^)  axes,  and  the  inclined  (x^ ' ,X£ ' ,x^ 1 )  axes.  These  entries  are 
defined  as  part  of  the  input  for  each  element.  Alternatively,  in  order  to  minimize 
the  amount  of  input,  the  option  has  been  provided  for  generating  internally  the 
direction  cosines  using  the  information  available  within  AD1NA.  Then,  denoting  by 
[e]  and  [ e ’ ]  the  strain  tensors  referred  to  the  global  and  inclined  axes, 
respectively,  the  transformation  of  strains  is  expressed  in  the  usual  form  as 


[e’]  =  [X]1  [e]  [X] 


BASIC  PROGRAM  MODIFICATIONS 


The  modifications  made  to  the  program  are  reflected  in  overlays  ADINA  and  ADINI 
and  to  a  greater  extent  in  subroutines  THDFE  and  SHELTH.  Minor  modifications  have 
also  been  made  in  some  other  subroutines. 

Modifications  made  to  ADINA  are  minimal  to  allow  for  the  addition  of  new  common 
blocks  associated  with  strain  computations  and  output,  and  to  perform  dimensioning 
and  initialization  of  certain  new  variables. 

Overlay  ADINI,  which  administers  input  and  solution  control  data,  was  modified 
by  allocating  permanent  memory  locations  at  the  beginning  of  the  blank  common  area 
for  variables  required  in  strain  computations  and  output.  Also,  certain  input 
statements  were  modified  to  accommodate  some  key  parameters  that  control  the  output 
of  stresses  and  strains.  In  particular,  two  parameters  were  added  to  the  input 
stream  to  serve  as  overall  control  flags  for  determining  whether  or  not  any  output  of 
stresses  and/or  strains  is  to  take  place.  These  represent  overall  control  parameters 
in  addition  to  the  ones  available  at  the  element  level. 


v/f"- ' 


l >  Tvi  Jii-.  i. 


•  *  m.  r 


±0.57738  02691  89626 


0.00000  00000  00000 
±0.77459  66692  41483 


Figure  2a  -  3  x  3  Mesh 


Figure  2b  -  4  x  4 
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Figure  2c  -  5  x  5  Mesh 


Figure  2  -  Standard  Meshes  for  Surface  Strain  Output 
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Subroutines  THDFE  and  SHELTH  were  modified  extensively  as  outlined  below: 

1.  New  common  blocks  were  added  to  allow  for  the  flow  of  information  between 
the  various  overlays.  Also,  dimension  statements  for  new  variables  were  incorporated, 
and  key  parameters  initialized. 

2.  Input  of  information  for  the  elements  was  expanded  to  include  certain 
control  parameters  that  monitor  the  input  and  output  functions  for  the  strains,  the 
direction  cosine  matrices  for  strain  transformations,  and  information  for  the  surface 
strain  output. 

3.  Because  the  strain  output  option  is  provided  for  a  limited  number  of 
elements,  arrays  were  generated  near  the  end  of  the  input  phase  in  both  subroutines 
to  store  the  above  information  and  other  pertinent  data  to  allow  for  the  subsequent 
retrieval  of  the  information  during  the  execution  (calculation)  stage.  Similarly, 
the  necessary • code  was  added  to  retrieve  the  information  needed  during  the 
calculations  of  stresses  and  strains  for  both  the  branches  dealing  with  the 
calculation  of  these  quantities  at  specified  nodal  points  and  the  branch  pertaining 
to  the  Gaussian  integration  points. 

4.  Additions  were  made  for  the  computation  of  the  coordinate  arrays  for  surface 
strain  output  locations,  particularly  for  the  case  when  these  locations  are  defined 
as  the  grid  points  of  standard  meshes. 

5.  Changes  were  incorporated  in  both  subroutines  to  compute  the.  strains  at 
points  other  than  originally  provided  for  in  ADINA. 

6.  The  computation  or  retrieval  of  the  direction  cosine  matrix  was  added  and 
the  transformation  of  strains  was  accomplished  based  on  standard  operations. 

7.  Several  control  statements  were  added  for  monitoring  output  at  the  element 
level  and  changes  were  made  to  accomplish  the  additional  output  functions  necessary 
for  the  transformed  strains  at  nodal  points  and  other  specified  locations. 

As  a  result  of  these  modifications,  the  following  output  options  can  be 
exercised  in  the  output  of  stress  and  strains. 

a.  The  output  of  stresses  and  strains  can  be  controlled  at  a  global  scale 
through  the  proper  definition  of  the  parameters  ITHSE  and  ITHSA,  respectively,  that 
allow  for  the  printout  or  suppression  of  the  output  of  these  quantities. 


b.  Control  of  stress  and  strain  output  at  the  element  level  is  accomplished 
through  the  proper  choice  of  the  value  of  the  control  variable  IPS,  originally  used 
to  exercise  the  option  of  printing  or  suppressing  the  output  of  stresses.  In  this 
respect,  the  option  of  outputting  stresses  and  strains  can  be  exercised  for  any 
number  of  elements  (belonging  to  various  groups  of  elements),  provided  these 
quantities  are  output  only  at  locations  specified  by  means  of  stress  tables  or  at 
Gaussian  integration  points.  The  strain  output  can  be  referred  to  any  arbitrarily 
oriented  orthogonal  coordinate  system  through  proper  input. 

c.  In  addition  to  the  above  option,  for  a  limited  number  of  elements  (that  may 
belong  to  different  groups  of  elements),  the  option  is  provided  for  the  output  of 
strains  at  specified  points  in  the  element  or  at  surface  locations  associated  with 
grid  points  of  3  x  3  to  5  x  5  meshes.  As  in  the  above  case,  these  strains  can  be 
referred  to  arbitrarily  oriented  coordinate  axes  through  the  input  of  direction 
cosine  matrices  or  to  the  natural  (r,s,t)  coordinate  system. 

EXAMPLE  PROBLEMS 

Some  simple  problems  were  analyzed  primarily  for  the  purpose  of  testing  the 
operation  of  the  modified  computer  program.  Brief  descriptions  of  the  example 
problems  analyzed  at  present  and  the  resulting  strain  output  are  given  below. 

PURE  SHEAR— NONLINEAR  STRAIN  CASE 

This  problem  represents  a  very  simple  example  of  a  nonlinear  strain  case  drawn 
out  of  Eringen.^  It  consists  of  a  cube  in  "simple  shear,"  as  shown  in  Figure  3a, 
where  Y  represents  the  angle  of  shear  and  SY  the  translational  displacement  in  the 
X-direction  at  a  distance  Y  measured  along  the  Y-r.xis,  so  that 


„  SY 
tan  =  —  = 


S 


(22) 


or 


Y  =  arctan  S 


(23) 
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The  Lagrangian  strain  tensor  E  representing  the  state  of  deformation  as  given  by 


Er ingen  is 


0  S/2  0 

S/2  S2/2  0 

0  0  0 


For  the  purpose  of  obtaining  numerical  results  using  the  modified  ADINA  program,  a 
single  three-dimensional  solid  element  with  21  nodes  was  used,  as  shown  in  Figure  3b. 
The  dimensions  of  the  cubical  element  were  chosen  to  be  1.0  in  all  three  coordinate 
directions.  The  structure  was  constrained  against  movement  in  the  Y-direction  so 
that  no  change  in  length  takes  place  along  the  Y-axis.  Also,  constraints  in  the 
Z-direction  were  imposed  and  the  movement  in  the  X-direction  of  nodes  2,  3,  4,  9,  10, 
11,  12,  17,  18,  19,  20,  and  21  (see  Figure  3b)  was  constrained  to  correspond  to  the 
movement  of  node  1.  At  node  1  a  prescribed  X-displacement  of  0.01  inch  was  imposed 
to  generate  the  desired  state  of  simple  shear.  A  nonlinear  static  analysis  was  then 
performed  using  a  material  model  corresponding  to  an  elastic-plastic  material  with 
isotropic  hardening  (Model  8)  with  the  constants  defined  as: 


Modulus  of  Elasticity  E  *  30  *  10 
Poisson's  Ratio,  V  =  0.30 

Yield  Stress  =  36  x  10' 

Hardening  Modulus  =0.0  psi 


psi 

psi 


A  step  size  of  0.1  was  used  leading  to  a  total  of  10  solution  steps.  The 
results  of  the  analysis  were  as  expected  so  that,  at  the  final  solution  step,  the  X- 
displacements  of  nodes  1  through  4  and  9  through  12  were  found  to  be  equal  to  0.01 
inches,  of  nodes  17  through  21  equal  to  0.005  inches,  and  the  remaining  displacements 
equal  to  zero.  With  reference  to  the  (r,s,t)  coordinate  system  shown  in  Figure  3b, 
the  diagonal  and  off-diagonal  components  of  the  Lagrangian  strain  components  were 
found  (at  the  end  of  the  10th  solution  step)  to  correspond  to  the  exact  values  of 


T  =  0.00005 

rr 

T  „  =  -0.01000 

rt 
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LINEAR  ANALYSIS  OF  A  STIFFENED  PLATE 


The  structure  considered  in  this  case  consists  of  a  cantilever  plate  with 

vertical  rectangular  stiffeners  along  its  longitudinal  edges,  as  shown  in  Figure  4a. 

The  plate  has  a  width  of  40  inches  and  a  thickness  of  1  inch.  The  stiffeners  are 

1  inch  by  6  inches  in  cross  section.  This  structure  was  previously  analysed  by 

2 

Bathe  and  Bolourchi  using  ALINA  SHELL  elements.  In  their  study  they  obtained 
favorable  comparison  with  the  solution  calculated  using  elementary  beam  theory.  In 
the  present  study,  the  same  structure  was  analyzed  using  three-dimensional  solid 
elements.  The  plate  and  stiffeners  were  each  modeled  using  two  three-dimensional 
elements,  for  a  total  of  six  elements,  66  node  points  and  108  degrees-of-f reedom. 

This  discrete  model  of  the  structure  is  shown  in  Figure  4b. 

A  linear  static  analysis  was  performed  considering  the  material  of  the 
structure  to  be  linear  isotropic  with  a  modulus  of  elasticity  of  30  x  10^  psi  and 
Poisson's  ratio  of  0.30.  The  two  concentrated  loads  acting  on  the  tip  of  the 
cantilever  were  taken  to  be  5,000  lb  each,  so  that  P  =  10,000  lb.  The  comparison 
of  the  vertical  deflections  of  the  tip  of  the  cantilever  obtained  by  means  of  beam 
theory  and  ADINA  is  given  as: 

Maximum  Deflection  -  Beam  Theory:  1.26  inch 
Maximum  Deflection  -  ALINA:  1.27  inch 

Typical  strain  outputs  for  elements  1  and  2  are  given  in  Table  1.  These  strains  have 

3 

been  multiplied  by  a  factor  of  10  in  order  to  retain  a  sufficient  number  of  signif¬ 
icant  figures  and  are  referred  to  the  (r,s,t)  local  coordinate  systems  shown  in 
Figure  4b.  The  strains  are  in  agreement  with  the  theoretical  values  based  on  the 
beam  theory  after  the  appropriate  transformation  of  the  strain  from  the  global  to  the 
local  system  of  axes,  as  verified  by  calculations  independent  of  the  computer 
program. 

CLAMPED  PLATE  SUBJECTED  TO  HYDRODYNAMIC  PRESSURE 

In  this  case  an  experimental  setup  of  a  structure  consisting  of  4  inch  steel 
plates  of  1  inch  thickness  and  clamped  at  both  ends  was  analyzed  using  the  modified 
ADINA  program. 
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E  : 

If  i 


Each  plate  basically  behaves  as  a  shallow  beam  of  4  Inch  width  and  1  Inch 
thickness  in  cross  section  fixed  at  both  ends  and  having  a  span  of  21  inches.  Due  to 
the  nature  of  tne  dynamic  loading,  this  problem  corresponds  to  a  nonlinear  fluid- 
structure  interaction  problem.  In  order  to  carryout  the  analysis,  the  modifications 
for  the  strain  output  were  incorporated  in  an  augmented  version  of  ADINA  that 
provides  the  capability  for  hydrodynamic  analysis  due  to  underwater  explosions. 

In  the  finite  element  model  of  the  structure,  a  total  of  20  three-dimensional 
solid  elements  were  used  leading  to  a  discrete  model  of  206  nodes  and  583  degrees- 
.  of-freedom.  A  linear  dynamic  analysis  was  performed  based  on  a  linear  isotropic 
material  with  the  modulus  of  elasticity  equal  to  30  *  10^  psl  and  Poisson’s  ratio 
of  0.30.  A  time  step  size  of  0.5  x  10”^  was  used  in  the  analysis  corresponding  to  a, 
total  of  100  time  steps.  The  time  history  of  strain  at  the  midpoint  on  the  top 
surface  of  the  beam  is  shown  in  Figure  5,  in  which  the  horizontal  scale  is  the 
longitudinal  strain. 
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Figure  5  -  Time  History  of  Strains  at  Midpoint 


CONCLUDING  REMARKS 

The  computation  and  output  of  strains  in  the  interior  and  on  the  surface  of 
finite  elements  is  of  significance  in  engineering  applications.  The  procedure 
outlined  in  this  paper  within  the  framework  of  ADINA  allows  for  the  computation  and 
output  of  such  strains  in  any  arbitrary  orientation.  Results  demonstrate  the 
applicability  and  versatility  of  the  method  developed  in  this  study  for  postprocessing 
of  the  strains. 
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